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We study the particle spectrum and the unitarity of the generic n-dimensional Weyl- 
invariant quadratic curvature gravity theories around their (anti-)de Sitter [(A)dS] and fiat 
vacua. Weyl symmetry is spontaneously broken in (A)dS and radiatively broken at the loop 
level in flat space. Save the three dimensional theory (which is the Weyl-invariant extension 
of the new massive gravity) , the graviton remains massless and the unitarity requires that the 
only viable Weyl-invariant quadratic theory is the Weyl-invariant extension of the Einstein- 
Gauss-Bonnet theory. The Weyl gauge field on the other hand becomes massive. Symmetry 
' breaking scale fixes all the dimensionful parameters in the theory. 

<n' 

<D 



INTRODUCTION 



General Relativity (GR) is expected to be modified (or perhaps replaced by a new theory) at 
both large and small scales. For large distances one can keep the theory intact with the high 
price of introducing a huge amount of dark matter and dark energy to explain the flattening of 
the galaxy rotation curves and the accelerated expansion of the universe. In the UV regime, even 
though there is no experimental necessity to modify GR, from the perturbative quantum field 
theory perspective, it is a non-renormalizable theory with no predictive power, hence it must be 
replaced by a UV-complete theory. 

There are many proposals to modify GR at both of the regimes, here we study the most general 
Weyl-invariant extension of GR augmented with the Weyl-invariant quadratic curvature terms built 
from the contractions of the Riemann tensor. Thus we both extend the diffeomorphism symmetry of 
GR to the local scale invariance and add quadratic terms which give better UV behavior. [Whether 
the theory really makes sense at the quantum level is an open issue and there is not much work 
in the literature on the quantization of Weyl-invariant gravity theories, see [l|], and the earlier 
references therein, for discussions of the one-loop beta functions of the Weyl-tensor square gravity 
and see 0] for a similar computation in the conformally coupled scalar-tensor theory.] Local scale 
invariance demands that the action has no dimensionful parameters, therefore Newton's constant 
^ 1 and any other dimensionful parameter, such as the mass of graviton, the cosmological constant 

appears only after the scale symmetry is broken (either by hand or spontaneously or radiatively). 
[Sec [3] for a recent extension of the Standard Model with local conformal symmetry.] For large 
distance modifications of GR, as an alternative to the dark energy (or cosmological constant), one 
could give a tiny mass to the graviton. But in generic n dimensions, there does not seem to be a 
consistent way to give a non-linear mass to the graviton within the frame of a single-metric theory. 
Ideally, one would like to find a Lagrangian which includes gravity and a Higgs type scalar field 
whose vacuum breaks the symmetry (not specified yet) and as a result graviton becomes massive. 
If one considers the rather special case of the n = 2 + 1 dimensions, such a theory exists: The 
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relevant symmetry is the Weyl symmetry [hence one not only has a scalar field but also a non- 
compact Abelian gauge field and a dynamical metric] and the action is given by the Weyl-invariant 
extension 0, @] of the new massive gravity (NMG) 0, 0]. NMG, currently, is the only known 
parity invariant, non- linear extension of the Fierz-Pauli massive gravity. This attractive feature of 
the model led to a recent activity in the 2 + 1 dimensional gravity theories 0-0. Weyl-invariant 
extension of NMG 0, 0| goes one step further and provides (in AdS not dS) an example of the 
mass generation for graviton via the symmetry breaking mechanism in perfect analogy with the 
Higgs mechanism in the Standard Model. For flat backgrounds, the symmetry (Weyl symmetry) is 
broken at the two-loop level [2^] through the Coleman- Weinberg mechanism 26] and the graviton 
again gets a mass. All the dimensionful scales in this theory are fixed by the symmetry breaking 
scale. 

In this paper we study the most general Weyl-invariant quadratic gravity theory in re dimensions 
with respect to its unitarity and stability and explore the perturbative spectrum about the (A)dS 
and flat backgrounds. Unitarity of the theory, that is the ghost and tachyon freedom highly 
restricts the spectrum and rules out a massive graviton for n > 3, therefore no viable extension 
of the Weyl-invariant NMG seems to exist beyond 2 + 1 dimensions. In fact, as we shall see, 
the only unitary quadratic theory in generic n dimensions is the Weyl-invariant extension of the 
Einstein-Gauss-Bonnet theory. 

The layout of the paper is as follows: In Section II, we briefly review the construction of the 
Weyl-invariant quadratic actions which was given in [J]. In Section III, we find the maximally 
symmetric vacua of the theory and expand the action around the vacua up to quadratic order in 
the fluctuations of the fields. Section IV contains the discussion of the decoupling of the fields 
and in Section V, we find the masses of the free fields and explore the unitarity regions. Some 
computations which we used in finding the second order expansion of the action are summarized 
in the Appendix. 



II. WEYL-INVARIANT QUADRATIC THEORY 

Generic re-dimensional Weyl-invariant quadratic theory was constructed in [J] which we briefly 



recapitulate here. [See also [27|, [28|] and see [29| for conformally invariant higher derivative scalar- 
tensor theories where the scalar field is coupled to the Euler densities without the Weyl gauge 
field.] The literature on the Weyl gauging (which was in fact the first example of upgrading a 



global symmetry to a local one) of the metric is vast, for a quick review see [3CJ, [3l|. To grasp the 
bare essentials of the Weyl gauging idea, consider the kinetic part of the scalar field action (with 
the mostly plus signature): 

s„ = -\f (PxV^gg^d^d^. (l) 



2 _ 

For n > 2 dimensions, this action can be made locally scale-invariant in a generally covariant way 
as follows 1 

9»u -> V = eK{x) 9^ * -> *' = e-*^*)*, (2) 

where C( x ) * s a rea l function and the partial derivative should be replaced by the gauge covariant 
derivative as 

7? — 2 

= - -z A^Q, (3) 



1 n — 2 case can also be handled similarly, but in this paper, we do not consider this rather unique case. 
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with A p being the non-compact Abelian gauge field transforming as 

A^A^A^-d^x). (4) 

Since we shall introduce dynamical gravity, we also need to know how the gauge covariant derivative 
acts on the metric. By simple inspection, one finds 

Dpgap = d^gap + 2A iX g a p. (5) 

By construction, these definitions yield 

= e 2C(x) 2W> (PM = e-^^V^, (6) 

where ([2]) and were employed. We also need to write a kinetic term for the gauge field which 
comes with a compensating scalar 

/, 2(n-4) 
d n x^g F^Fi", (7) 

where F^ v = d^A u — d^A^ is the usual gauge invariant field strength, and a dimensionless param- 
eter e is introduced (unitarity requirement will later restrict it). Note also that for n = 4, the 
compensating factor drops out since the Maxwell action is already conformally invariant. Weyl 
invariance in the gravity part follows with the help of a Weyl- invariant "Christoffel connection" 
defined as 

rj„ = \g Xa (?W + v»9iut ~ v °9»v) ■ (8) 

From this connection a Weyl-invariant Riemann tensor can be defined as 

#Vbi A \ = W CT - d a v % + r^f I. - r^r^ p , (9) 

which, after making use of ([5]) becomes 

where we have used the notation 2Ar p B a i = A p B a — A a B p . Here is the usual covariant derivative 
compatible with the Riemannian metric g pv : For example one has V p,A v = d^Ay — T pu A a . Note 
that pU|) does not have the original symmetries of the Riemann tensor but this is not needed, all 
we want is to find a short way to construct Weyl-invariant actions for generic higher curvature 
gravity theories. For this purpose (|l(jp can be used. After contraction, Weyl-invariant Ricci tensor 
can be found as 

R„ a [g,A] = Rf* vtl(r \g,A] 

r 9 i ( n ) 

= R vtT + F va - (n - 2) [V a A u - A v A a + A l g va \ - g ua V ■ A, 

where V • A = V ^A^ . One more contraction gives the scalar curvature 

R[g, A] = R — 2(n — 1)V ■ A — (n— l)(n - 2)A 2 , (12) 

which transforms like the inverse metric as (R[g,A])' = e~ 2 ^R[g, A]. Weyl invariance allows a 
potential term for the scalar field yielding the action 

S * = ~\f d n x^{v^V^ + (13) 
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where v > is a dimensionless coupling constant. Positivity of v is required for the existence of 
a ground state, further constraints can come from the requirement of unitarity or the existence of 
a maximally symmetric vacuum, as we shall discuss later. Collecting all the pieces together, the 
generic Weyl-invariant quadratic action is given as 0] 

Swi = J d n x^g~ ^ a $ 2 R + <S>^r [aR 2 + 0^ + j + 5$ + S A » , (14) 

where Sap and S$ are the Weyl-invariant actions for gauge field and the scalar parts given by 
(J7]) and (fT3j) . The action (fT4"|) is invariant under the transformations ([5]) and d3]). Note that Weyl 
invariance alone leaves 7 free dimensionless parameters, one of which can be eliminated by scaling 
the total action. Here we have done the scaling in such a way that the kinetic part of the scalar 
field action comes with the canonical 1/2 factor. 

Apparent simplicity of (|14|) is somewhat deceptive since the explicit forms of the curvature 
terms yield a highly complicated action in terms of the fields [g^jA^, $]. The explicit form of the 
scalar curvature square follows from (|12|) 

R 2 =R 2 - 4(n - l)i?(V • A) - 2(n - l)(n - 2)RA 2 

+ 4(n - 1) 2 (V • A) 2 + 4(n - l) 2 (n - 2),4 2 (V • A) (15) 
+ (n-l) 2 (n-2) 2 >l 4 , 

where V.A = V M ^, ^4 2 = A^ and A 4 = A^A U A V . The square of the Ricci tensor follows 
from (HI|) 

^jL = ~ 2 ( n ~ 2)1^ VM** - 2R(V ■ A) + 2(n - 2)R ftu A^A v - 2(n - 2)iM 2 

- 2(n - 2)F^V^ + F 2 , + (n - 2) 2 (V^ M ) 2 + (3n - 4)(V • A) 2 - 2(n - 2)%^^ 
+ (4n - 6)(n - 2)A 2 (V ■ A) + (n - 2) 2 (n - 1){A) A . 

(16) 



Finally, the square of the Riemann tensor follows from (|1U|) 

2jLpr = " 8fl^ V A" + 8/^ AM" - ARA 2 + nF 2 v -, , w , _ M , ,, 



fliUr = - 8^V"A" + 8/^AM" - ARA 2 + nF 2 „ + 4(n - 2)(V M A,) 2 + 4(V • A) 2 

+ 8(n - 2)(A) 2 (V • A) - 8(n - 2)^^VM" + 2(n - l)(n - 2)(^) 4 . 



(17) 



Inserting (|15p , ()16|) and (jl7j) into (|14|) yields an action in the Jordan frame with many non-minimal 
interaction terms: For example the gauge field is allowed to self interact at the quadric A 4 level, 
in a gauge invariant way, even though this is an Abelian gauge theory. 



III. THE ACTION UP TO QUADRATIC ORDER IN THE FLUCTUATIONS OF THE 

FIELDS 

Now that we have constructed the most general Weyl-invariant quadratic gravity theory, we 
can proceed to find its perturbative spectrum, namely the particle content of the theory about 
its maximally symmetric vacuum or vacua. Naively, one should find the field equations first and 
find the constant curvature vacuum solutions, then linearize the equations about any one of these 
vacua and identify the propagating degrees of freedom. But, this is a very complicated procedure, 
even the field equations are hard to find (see 0] for the field equations in the n = 3 case). Instead 
of this cumbersome procedure, we will follow the technique used in [33| and directly expand the 
action about its assumed maximally symmetric constant curvature vacuum up to second order in 
the fluctuations of the fields. This method not only will give the relation between the parameters 
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of the theory in the vacuum, but it also will simplify the computation of the parts quadratic in the 

I L 

fluctuations of the fields. For n = 3, we have given the computation in [5]. Note that one does not 
need to go to the Einstein frame to determine the particle spectrum of the theory, one can directly 
work in the Jordan frame. Equivalence of the specific case of the conformally coupled scalar tensor 
theory with regard to its particle spectrum was shown in [jj. 

Let us consider a dS or an AdS background (flat vacuum can be obtained in the limit) for which 
the vacuum values of the fields are 

<*W = m (n " 2)/2 , A» ac = 0, g^ = 9»u, (18) 

where m is of the mass dimension and appears either by the requirement that an (A)dS vacuum 
exists or as we shall briefly discuss later, in the flat space case, conformal symmetry is broken 
at the loop level via the Coleman-Weinberg mechanism 0,0 where the scalar field acquires a 
non-zero vacuum expectation value. Hence in either case symmetry is broken in the vacuum (|18p . 
[There is a caveat here: Strictly speaking, the Coleman-Weinberg potential was computed in three 
[25| and four dimensions [2^], for larger dimensions we expect the symmetry to be broken at the 
loop level, but no explicit computation is available.] To study the particle spectrum of the model 
we need to consider fluctuations about the vacuum as 

$ = m ("- 2 )/ 2 + T $ L , A lt = TA%, g^ = g^ + rh^, (19) 

where we have introduced r, a small dimensionless parameter to keep the track of the expansion 



orders, at the end it can be set to unity. In what follows, we will use the conventions given in 331 ] . 
Let us represent the expansion of the action (|14p under the fluctuations (|19p as follows: 

S WI = J <Px^}\c{t°) + rCir 1 ) + r 2 £(r 2 ) +•••}, (20) 

where C(t°) is a constant built from the vacuum values of the fields which is irrelevant for the 
perturbative spectrum that we shall discuss. On the other hand C{t 1 ) reads 

£( r !) = (^- m ^<$> L + ^m n - 4 h) (CA 2 + 4a Am 2 - urn 4 ) , (21) 

where 

Criticality of the action for arbitrary variations of the fields about the vacuum, from (|2ip . yields 

CA 2 + 4oAm 2 - vra 4 = 0, (23) 

which reduces to the one given in 0, IB] for n = 3 and for the NMG condition (7 = and 
8a + 3/3 = 0). One can read (|23p in two different ways: First, one can assume that the symmetry is 
broken (namely m is given), then A is fixed by the symmetry breaking scale. Or, one can assume 
that A is given (meaning that the existence of (A)dS vacuum is imposed) and m is determined. 
Let us discuss the first case (the second case follows exactly the discussion in 0]). For n / 4, there 
are two vacua given as 



2m 

A± = 



C 



Cu 

T 



n ^ 4. (24) 
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As long as a 2 + ^ > 0, there is at least one constant curvature vacuum. Four dimensions is unique 
in the sense that there is a single vacuum 



A 



vm 



(25) 



Finally, we need to find the C(t 2 ) part of the action, which requires a rather long computation. 
A naive expansion of the quadratic parts of the action would result in a cumbersome expression 
which will not be explicitly background diffeomorphism invariant. To simplify the computation, 
one should use background diffeomorphism invariant quantities, such as the linearized Ricci scalar 
or the Einstein tensor as well as the self-adjointness of the involved operators and the linearized 
Bianchi identities. We give some essential parts of this computation in the Appendix and write 
the result here 



£(r 2 ) 



+ m 2 



1 



-m 



A 



n—iy^av 



2A 



An 4 n 
-a + r/3 



-7)AS£, + (2a + + 2 7 )(^D - V u V u )r l 



n-2 



C^ + 2a 



n — 2 n — 1 n 

2 f„^n\ , i 



2a 



n — 1 



n 



2 n .L 
jiv 



Rl$l + 



n [n(n-6)C K l . (n + 2) 2 
+ 4crA - ^-m u 



2(n - 2) L (n-2) m 2 



n-2 



$2 _ (5 $ L )2 



m 



„n— 4 



4(n - l)a + n/3 + 47 



+ m 



n— 4 



1 



A 



2(n-l)C^ + 4<r(n-l) + 



2 v 
n- 2- 



4(n - l) 2 a + n/3 + 4 7 J (V • A L f + -m n ~ 4 [(n 2 - 2n + 2)(3 + 2(3n - 4) 7 + 2e (F ( 



L \2 



2m 



n-2 



(2n[n - l)a + (3n - 4)/3 + 87) — =■ + 



(n-l)(n-2) (n-2) 
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4 2 



where fl^ and Q^ v are the linearized Ricci scalar and Einstein tensors defined by 34], \3 



(26) 



1 



-,9nvR 



n-2 ' 
2A 
~ n-2 



(27) 



ah 



[IV 



In ([2"Bj) . vacuum equation ((23D has been used. But, (f2"B]) is still a very complicated coupled system 
of relativistic "harmonic oscillators" which need to be decoupled before one can discuss what the 
true propagating degrees of freedom are. It pays to fix the Weyl gauge which we shall do first 
below. 



IV. DECOUPLING OF THE FIELDS 
A. Gauge fixing condition 

To remove the redundancy we need a proper gauge-fixing condition. Let the gauge-covariant 
derivative act on the gauge field as: 



V a A v = V n A u + (n - 2)A n A v . 



(28) 
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Note that, compared to (0), we used V M in the first term instead of the partial derivative. Under 
the transformations ([21 ED , it is easy to show that the divergence transforms as 

(V^)' = e~ 2 <W {v^ _ V^C{x)) . (29) 
Therefore, we can choose a Lorenz-like condition 

V^A* = V • A + (n - 2) A 2 = 0. (30) 

It is important to note that V^d^C, = is also Weyl- invariant. [This is a Weyl-invariant generaliza- 
tion of the leftover gauge-invariance, d 2 C = 0, after the usual Lorenz gauge d^A^ = is chosen.] At 
the linear level, (|3U|) reduces to the background covariant Lorenz gauge fixing condition : V-Al = 0. 
With this choice, three terms in the action (|26p drop out. 



B. Redefinition of the tensor field 



In order to decouple the tensor and scalar fields, let us choose 

~ 4 2-n 

/V = V' ~ n -2 m 2 9 ^ Ll ( 31 ) 

which yields the following relation between the linearized curvature expressions which we all need 
in the computation 

R%, =R% + ^^ m ^ {( n ~ 2 )v,A$l + 9^ L ) , 

~ 4 2-n / , , 2n \ 

Rl =Rl + n - + -A$ L , 

n — 2 V n — 2 > 

Q% =9%, + 2m 2 ^r(v fl d u <S> L - g^U^ L - 2(n - 2)Ag^$ L ), 

h ,UgL u + 4m 2-n [ m ^R^ L + 2 ^l$ iD $ L + J^^^l) , 

K»? =(2iv? + 4m 2 -((n - 1)(D$ L ) 2 + _^L_ A 2 $| + ^^Wn^) 

+ 2m ^ ((n - 2) J R L D$ L + 2A J R L $ L ) . 

After making use of (|31|) and the gauge-fixing condition and the vacuum equation (|23p (which 
removes the $\ term), the vector field is decoupled from the rest and the quadratic part of (fT4"|) 
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boils down to a more transparent form 



Swi 



d n x, 



,n— 4 



g< m 



+ [a 
1 r 16 



/ 2nA 2A 

n« + nP 



>n - 2 
n - 4 



n 



4(n-4)A ,rn 2 
(n-l)(n-2f + ^ a)h 



L 



(n - 2)2 



+ 



16 (n 



m 2 (n — 2 
6 n 



+ 8m 2 



1 

H — m 
2 



n— 4 



2m 



n-2 



"(n 2 - 
(2n( 



4 

^2n(n 

a + 
a + 



/3-(n-3) 7 )i?! + (/3 + 4 7 )(^) 2 
- l)a + (3n - 4)/3 + 87) A + 

(n 



-<7 + l 



m z 



2)2 l 



n 



4(n- 1, 
n 

4(n- 1) 



1 



P + 



n — 1 
1 



n 



1 



L \2 



2n + 2)/3 + 2(3n - 4) 7 + 2e (F^) 



l)a + (3n - 4)/3 + 87) A + 



(n-l)(n-2) (n-2) 



16 



(33) 

We still have the coupling between the tensor and the scalar fields. But this is not a problem 

since it drops out once unitarity is imposed to remove the higher-derivative Pais-Uhlenbeck term 
(□$£) 2 . It is remarkable that both of these unwanted terms come with the same coefficient, which 
must be set to zero: 



a + 



4(n - 1) 



1 



11 



-7 = 0. 



(34) 



It is worth mentioning that this condition is a necessary but not a sufficient condition for the 
unitarity of the theory. In three dimensions this condition gives the NMG theory after noting that 
the Riemann tensor can be written as R 2 avpa = 4i? 2 — R 2 , which is the Gauss-Bonnet identity, 
therefore it is easy to see that (f3"4"|) reduces to 8a + 3/3 = once the square of the Riemann tensor 
is eliminated with the help of the Gauss-Bonnet identity. Now that all the fields are decoupled, we 
can find the masses and discuss the unitarity regions in the theory. 



V. PARTICLE SPECTRUM AND THEIR MASSES 



In order to read the masses of the fields, let us write fl3"3"D as 



where 



n—A 



Swi 



2nA 2A n 
a + -/3 



4(n - 4)A 



+ Q 



>n-2 n-2 
n — 4 



/3-(n-3) 7 m + (P + 4ry){Gt,) 



n—A 



(n 2 - 2n + 2)f5 + 2(3n - 4)7 + 2e| (FM 2 



2m 



n-2 



f 2n(n - l)a + (3n - 4)/3 + 87) —= + 



A (n-l)(n-2) (n-2) 



(7 + 



16 



A 2 



16 



(n - 2)' 



■ (2n{n - l)a + (3n - 4)/3 + 87) 4r + 8^— i(<7 + ll (d M $ L ) 
V / m z n — 2) J 



(n-1) 



(n-2) 



(35) 



(36) 
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The gauge field is of the massive Proca type and the scalar field just has the kinetic term. On 
the other hand, the tensor part is still complicated. Let us first look at the scalar and gauge field 
parts. 

One can choose e in such a way that the kinetic term of the gauge field becomes zero which 
renders the gauge field non-dynamical, but instead we choose 

e = -i((n 2 -2n + 2)/3 + 2(3n-4) 7 + l/2), (37) 

1 . (n-4) 

to normalize the kinetic part to its canonical value —4 after scaling the field by a factor of m 2 . 
We thus have 

Ca, = -\(F^?-\m 2 a AI, (38) 

where the mass-square reads 

r 1 r (n — 2) 2 i 

M\ =4(n-4)[2(n- l)o + /3jA+ \ 2{n - l){n - 2)a + y - — ^^\m 2 . (39) 

In getting this equation, we have imposed the constraint (|34|) which was necessary for unitarity. 
Note that, for the gauge field alone M\ > is sufficient for unitarity and in this case the theory 
has a massive (or massless if the bound is saturated) spin-1 excitation. But, below we will see that 
there will be another relation between a, f3 and 7 coming from the unitarity of the spin-2 part. In 
three dimensions, (|5U|) reduces to the one found in 0]. 
On the other hand the scalar field part becomes 

AM 2 1 

c * = - (^)W 8 *^ 2 - (40) 

which does not impose any new condition: The condition that makes the gauge field non-tachyonic 
makes the scalar field non-ghost. Also, when the gauge field has zero mass, the scalar field becomes 
non-dynamical. For other values of M|, the scalar field can be re-scaled to have = — ^(c^^l) 2 . 

Finally, let us investigate the spin-2 part. The Lagrangian density, following the procedure of 
0, can be written in terms of two auxiliary fields <p and f^ u , as 

/V =ohrg L ^h) + bR\ + c(Q L ^(h)f 

2 2 (41) 

^ahrgfah) + roiM + vRl - - - n 

where / = g^ u and 

„_,/ 2nA 2A n 4(n-4)A m 2 

a 55 ~ m \^2 a + ^ ~ (n - l)(n - 2) 7 + T* 

n-4 \ (42) 



4 

c=m n - 4 (/3 + 4 7 ). 

We have to first determine mf in terms of the parameters of the theory. This can be done by using 
the field equations for the auxiliary fields 

/p^W+ffp^fii. <P= ^Rl, (43) 



10 



in (|41|) to obtain 



•L ,,> , 1 ,„£ „„2 / 1 (» - 2) 2 N o2 



which yields 



1 Ah ^ Jn^2) 2 

— 7 , and o = T) — — 

m,2 2m( 4(n — l)m^ 



2' and 6 =TT-2-77Z 7TZ2" ( 45 ) 



Then ([41 [) reads as 

c K „ = ( .*r + os» + ^ - 46+c(n _ t ' 2 2 )V(n _ 1) - - n («) 



The unitarity condition (|34p in this parametrization becomes 

46 + c^-^ = 0, (47) 
n—1 

which decouples (p and drops the term <pRl in (|46|) . Therefore, the gravity part becomes: 



l Kv = {ahr + DQiM - ^(/J, - / 2 ). (48) 



In order to obtain the spectrum and read the masses, one should decouple the h^ u and fields 
which can be done with the following redefinition (assuming a^O, see below for the a = case) 

1 

Inserting this in ([48 p results in: 



4w = «^^(h) - ^ro^U) - ^(J% - f 2 )- ( 5 °) 



The first term is just like the linearized part of the pure (cosmological) Einstein-Hilbert theory 
with an effective Newton's constant, therefore it propagates a massless spin-2 field as long as a < 0. 
The second part is the Lagrangian density of a massive Fierz-Pauli spin-2 field. Because of the 
coefficient of the kinetic part, it is clear that massless and massive spin-2 fields cannot be unitary 
at the same time. The only solution is to freeze the massive spin-2 field by giving it an infinite 
mass: c = 0, which, together with our earlier unitarity condition (|34[) . yield 

4 7 + p = 0, a = 7. (51) 

This is exa ctly the Einstein-Gauss-Bonnet theory which has a massless unitary excitation in (A)dS 
as long as [32| 

4(n-3)(n-4) 7 A . 
a> (n-l)(n-2)m2 ' 

Therefore in n > 4 dimensions, out of all the Weyl-invariant quadratic gravity theories, unitarity 
condition singled out the Weyl-invariant extension of the Einstein-Gauss-Bonnet theory. The mass 
of the gauge field in this theory reduces to 



M 



8(n - 3)(n - 4) 7 A + ft 1 ^ + 2(n - l)(n - 2)a)m 2 . (53) 



11 



Note that M\ > 0, a condition for the unitarity of the scalar and gauge field parts, gives 

4(n-3)(n-4) 7 A _ n-2 
~ (n- l)(ra-2)m 2 8(n-l)' 1 j 

which is a weaker condition than (|52|) . Finally one can check that the stronger unitarity condition 
(|52|) . which now becomes <r > — is compatible with the existence of a maximally symmetric 
vacuum ([25]) for both dS and AdS spaces. 

Consider now the critical case, a = in (1481) for which one has 



^ = ^i(/)-^(/i-A (55) 

where we have used the self-adjointness of the involved operators in the first term. Now variation 
with respect to h^ u yields 

GiM) = 0, (56) 

which can be solved as 0] 

f l u, = V lt B v + V v B lt . (57) 
Inserting this into (|55|) gives up to boundary terms 

where F^u is the field strength of the S-field. For c > 0, ([58]) describes a massive spin-1 excitation 
with M 2 = in analogy with the critical point of NMG 0]. One can check that the unitarity 
conditions found before are compatible with the criticality condition and the full Weyl-invariant 
quadratic theory propagates a unitary massive Weyl gauge field and a massless scalar field in 
addition to the just discussed massive spin-1 field in dS. 



A. Flat background for n = 4 

In the flat backgrounds in four dimensions, as we noted before, the symmetry is broken via 



the Coleman- Weinberg mechanism 26] for which the one-loop effective potential becomes V($) = 
ci<I? 4 (log(3>/m) + C2), where the actual values of the constants are not relevant here: All that 
matters is that the scalar field gets a vacuum expectation value and the symmetry gets broken. 
Then ([M]) together with (3 + 47 = gives 

?7l 2 

^h^ = - ~y^ v q^ u , 

^=-\{B%f -i(l + 12a)m 2 A|, (59) 

C^ = -^(l + l2a){d^ L f. 

So in the Weyl-invariant extension of the Einstein-Gauss-Bonnet theory in four dimensions there 
is a unitary massless spin-2 field, a massive spin-1 and a massless spin-0 field for a > 0. Newton's 
constant and the mass of the gauge field are fixed by the vacuum expectation value of the scalar 
field. [Note that any other n > 4 theory has the same spectrum.] 
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VI. CONCLUSIONS 

We have found the perturbative particle spectrum and discussed the unitarity of the n- 
dimensional Weyl-invariant quadratic gravity, constructed in [ij], about their (A)dS and flat back- 
grounds. Three dimensional case, that is the Weyl-invariant extension of the new massive gravity, 
studied at length in [5], has the remarkable property that through the breaking of the symmetry, 
the graviton gets a unitary Fierz-Pauli type mass. The hope was to extend this mass-generation 
mechanism to four and more dimensions, but, as we have shown in this paper, the only unitary the- 
ory, beyond three dimensions, among the Weyl-invariant quadratic theories is the Weyl-invariant 
Einstein-Gauss-Bonnet model which propagates a massless spin-2 particle as well as massive spin-1 
and massless spin-0 particles. [ For an interesting Higgs-type mechanism in the first order formal- 
ism of gravity see [361 ] .] Now that we have shown the unitarity of the theory, it would be interesting 
to study the Weyl-invariant Einstein-Gauss-Bonnet theory in four dimensions with respect to its 
black hole and cosmological solutions. 



VII. APPENDIX 



We would like to expound upon our computation of the second order action in the fluctuations 
of the fields here. The expressions below are valid up to 0{t 2 ). The scalar field action reads 



5$ = -~ j d n x^f^g\vm n + r 



m . n+2 n 

h + 2m 2 <£ L 

2 n-2 



{d^L? + (n - 2)m— V • A L <5> L + 



(n-2) 2 _ 2 , 2 



-m n -'Aj 



(60) 



n(n + 2) n _ 9 



71 n + 2 



m 



m 



n _ 2 m 2 t ,h* L + —vh i -—vh* v 



(n - 2) 2 

The coupling term between the scalar field and the Ricci scalar yields up to quadratic order 
J d n x^$> 2 R = 
d n x^m n - 2 - 



2n 



n-2 
1 



A + r 



n a , ^ 4n . 2^ 

-Ah + R L ^ -Am 2 

n — 2 n — 2 



1 



— A/i 2 „„ - -Ah 2 - hR L - J^Am 2 -"$ L 
n-2 ^ 2 n-2 



(61) 

To expand the quadratic curvature parts of the action, the following steps proved useful: First let 
us write 



I 2 



d n x^g§^=^ [aR 2 + (3Rl u + 7 # 2 
cPx^gi (m*? + r<D L ) ^ (X + rl« + r 2 X^) }, 



(62) 



where one finds 



X 



nC 



A 2 



'2(n-4) ' 



nC 
: 4(n-4) 



A 2 h + 



(n - 2)C 
2{n - 4) 



(63) 



ARj 
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Now in order to find X^ 2 \ we rewrite it as: 



2 

fiupa 



(2) 



-g ((a - 7 )i? 2 + (/3 + 4 7 )i?^ + 7XB 



(2) 



(64) 



where = Rp\upa ~ 4-R 2 j, + -R 2 is the Gauss-Bonnet combination. From [34|, 135J], one can write 
the as: 



2 



4nA 4A „ 
■a + T f3 



n 



n — 1 



8A 
n — 1 



+ (2a + /3 + 2 7 ) □ - V M V,)i? L + (2a + — ^-/3 - 2(n % )g R L 

- ' n-z\ n — 2 n — 1 / 

(65) 



+ (/3 + 4 7 ) + ^A 2 V " ^A 2 ^/i 
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